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We point out that the decay modes D^Dj and D^D^ have no 

factorized contribution. At quark level these dacays can only proceed through the 
annihilation mechanism, which in the factorized limit give zero amplitude due to 
current conservation. In this paper, we identify the dominating non-factorizable 
(colour suppressed) contributions in terms of two chiral loop contributions and one 

> : 

CN . soft gluon emission contribution. The latter contribution can be calculated in terms 

^ , of the (lowest dimension) gluon condensate within a recently developed heavy-light 

' chiral quark model. We find bracking ratios BR(B^ — > DfDj) ~ 7 x 10^^ and 



BR{B'^ D+D-) ~ 1 X 10"^ 



INTRODUCTION 



X 

^-^ ' There is presently great interest in decays of i?-mesons, due to numerous experimental 



results coming from Ba Bar and Belle, and later at LHC 

It has bee„ shown Q that so.e classes of «-,„e.o„ decay amplitude. exKib. QCU fac 

torization. This means that, up to a^/vr (calculable) and AqcD/fnb (not calculable), their 
amplitudes factorize into the product of two matrix elements of weak currents. Typically, 
the decay amplitudes which factorize in this sense are B ^ and B Kir where the 
energy release is big compared to the light meson masses. However, for various decays of 
the type B — > DD where the energy release is of order 1 GeV, QCD factorization is not 
expected to hold. (Here B, D, and D contain a heavy b, c, and anti-c quark respectively). 
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|aksel.hiorth@fys.uio.no| 



2 



Such decay modes have been considered in connection with intermediate DD states for other 
i?-decay modes P]- 

In a previous paper ^j, it was pointed out that the decay mode K^K^ was zero 

in the factorized hmit due to current conservation. However, there are in that case non- 
factorizable (colour suppressed) contributions in terms of chiral loops and soft gluon emission 
modelled by a gluon condensate. 

In this paper we report on the following observation: The decay modes B'^ D^Dj 
and — > D^D~ have no factorized (colour non-suppressed) contributions. At quark level, 
these decays a priori proceed through the annihilation mechanism hs —>■ cc and bd — > cc, 
respectively. However, within the factorized limit the annihilation mechanism will give 
a zero amplitude due to current conservation, as for K^K^. But there are non- 

zero factorized contributions through the axial part of the weak current if at least one 
of Z^-mesons in the final state is a vector meson D*. Such contributions are, however, 
proportional to the numerically non-favourable Wilson coefficient Ci, which we will neglect 
in this short paper. In contrast, the typical factorized decay modes which proceed through 
the spectator mechanism, say B'^ D^Dj , are proportional to the numerically favourable 
Wilson coefficient C2. If the mesons in this amplitude are also allowed to be vector mesons, 
such amplitudes will generate non-factorizable (~ 1/Nc) chiral loop contributions to the 
process B^ D+Dj^ due to /('''-exchange. These will be considered in the present paper. 

There are also non-factorizable (~ l/iVc) contributions due to soft gluon emission. Such 
contributions can be calculated in terms of the (lowest dimension) gluon condensate within 
a recently developed Heavy Light Chiral Quark Model (HLxQM) 4], which is based on 
Heavy Quark Effective Theory (HQEFT) [5]. This model has been applied to processes 
with 5-mesons in The gluon condensate contributions is also proportional to the 

favourable Wilson coefficient C2. 

In the next section (II), we shortly present the four quark Lagrangian at quark level. 
In section HI we present our analysis of chiral loop contributions within the heavy light 
chiral perturbation theory. In section IV we give the calculation of non-factorizable matrix 
elements due to soft gluons expressed through the (model dependent) quark condensate. In 
section V we give the results and conclusion. Throughout the paper, we will give formulae 
and figures for the decay mode B^ DfDj. The treatment of 5° D^D^ will proceed 
analogously. 
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II. EFFECTIVE NON-LEPTONIC LAGRANGIAN AT QUARK LEVEL 

Based on the electroweak and quantum chromodynamical interactions, one constructs an 
effective Lagrangian at quark level in the standard way: 

Cw = Y.^^{^i) Q^{^i) ^ (1) 
i 

where all information of the short distance (SD) loop effects above a renormalization scale 
/i is contained in the (Wilson) coefficients Cj. In our case there are two relevant operators 

Qi = 4(g^7"6i) {clIcCl) , Q2 = 4 (c^t^^l) (^lT^cl) , (2) 

for q = d,s. Penguin operators may also contribute, but have small Wilson coefficients. We 
may write Cj = —^VcbV*qai , where ai ~ 10~^ and 02 ~ 1 at the scale = rrib. Per 
ing perturbative QCD corrections within Heavy Quark Effective Theory (HQEFT) 
effective Lagrangian (Q) can be evolved down to the scale fi ^ ~1 GeV jsll^, where one 
finds |ai| ^ 0.4 and \a2\ — 1.4. The b, c, and c quarks are then treated within HQEFT. 

In order to study non-factorizable contributions at quark level, we may use the following 
relation between the generators of SU{3)c I, n are colour indices running from 1 to 3): 

SijSln = J^Sin^lj + 2 , (3) 

where a is the color octet index. Then the operators Qi^2 Kiay, by means of a Fierz trans- 
formation, be written in the following way : 

Qi = + 2Q2 , Q2 = ^Qi + 2Q, , (4) 

where the operators with the "tilde" contain colour matrices: 

Qi = 4(g^7°t'^6i) {clIJ^Cl) , = 4 {cL^t^L) {qLlaf'cL) . (5) 

To obtain a physical amplitude, one has to calculate the hadronic matrix elements of the 
quark operators Qi within some framework describing long distance (LD) effects. 

As an example of a typical factorized case we choose the amplitude for 
obtained from ((H) and (0): 

{D:D+\CwW)f = -{C2 + ^C,){D:\s^^^,c\0){D+\c^^b\W) , (6) 



FIG. 1: Factorized contribution for D^Dj through the spectator mechanism, which does 

not exist for decay mode DfD~ we consider in this paper. The double dashed hnes represent 

heavy mesons, the double lines represent heavy quarks, and the single lines light quarks. 
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FIG. 2: Factorized contribution for B^ DfDg through the annihilation mechanism, which give 
zero contributions if both and D~ are pseudoscalars. 

which will in section III be compared with our chiral loop contributions. This term is 
proportional to the D-meson decay constant times the Isgur-Wise function (for B ^ D 
transition) and is vizualized in figure ^ 

The factorized amplitude for 5° DfD~ obtained from ((1} and is vizualized in 
figure El and is given by 

{D;Dt\CwmF = 4(Ci + ^C2){D;Di\T^^,CL\0){0\TLYbL\B^) . (7) 

Unless one or both of the D-mesons in the final state are vector mesons, this matrix element 
is zero due to current conservation: 

{DtD;\c^^c\0){0\dri,b\W) ~ fsiPD+PDr {DtD;\c^^c\0) = . (8) 

The genuine non-factorizable part for B'^ DfD~ can be written in terms of coulored 
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currents (see eqs. (jH) and (0)): 

{D;Di\CwmNP = 8C2 {D;Dt\iqLrt%L) {cLlaf'cL) m (9) 

We observe that the annihilation mechanism amphtude in the non-factorizable case has 
the numerically favourable Wilson coefficient C2. This amplitude is, within the HLxQM 
vizualized later in figure |3] . 



III. HEAVY LIGHT CHIRAL PERTURBATION THEORY 

Our calculations will be based on HQEFT [s^, which is a systematic l/mq expansion 
in the heavy quark mass ttlq. The heavy quark field Q{x) = b{x) (eventually c{x) or c) is 
replaced with a "reduced" field Q^\x) for a heavy quark, and Qi^\x) for a heavy antiquark. 
These are related to the full field Q{x) in the following way: 

Q(±)(a;) = P±e^^-«''-Q(x), (10) 

where P± are projecting operators P± = (1 ± 7 ■ t')/2. The Lagrangian for heavy quarks is: 



Chqeft = ±Q1^^ tv ■ D Q(±) + OimQ') , (11) 

where is the covariant derivative containing the gluon field. In [3| the I/tuq corrections 
were calculated ioi B — B -mixing. In this paper these will not be considered. 

Integrating out the heavy and light quarks, the effective Lagrangian up to 0{'m^) can 
be written as 0, [l^ 



H'it\v ■ V,^'^^^ - g^Tr [Hj,"^^ hI^^^^^^.A^ , (12) 
where Hi"^^ is the heavy meson field containing a spin zero and spin one boson: 

Hi^^ ^P±iPi^h'-^Pl?l.) . (13) 

The fields pI^\p1^j ^) annihilates (creates) a heavy meson (vector for M = /i and pseu- 
doscalar for M = 5) containing a heavy quark (anti-quark) with velocity v. Furthermore, 
^■^ba ~ i^baD'^ ~ cind o, h are fiavour indices. The vector and axial vector fields and 
contain the field ^ (and its first derivative) which is a 3 by 3 matrix containing the (would 
be) Goldstone octet (vr, rf) : 

V,^\{ed,i + id,e) ■ A,^ -'-{i^d,i~id,e) , i^expi{Jl/f) (14) 
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where / is the bare pion couphng, and 11 is a 3 by 3 matrix which contains the Goldstone 
bosons 7T,K,ri in the standard way. The axial chiral couphng is Qj^ ^ 0.6. Eqs. ()12|). (|T^. 
and (fT^ will be used for the chiral loop contributions. 

The simplest way to calculate the matrix element of four quark operators like (^1,2 in 
eq. (1) is by inserting vacuum states between the two currents, as indicated in section II. 
This vacuum insertion approach (VSA) corresponds to bosonizing the two currents in Q12 
separately and multiply them, i.e. the factorized case. Based on the symmetry of HQEFT, 
the bosonized current for decay of the bq system is jj, [l^ : 



WYQi"^ — ^Tr[eV^i^r^J , (15) 

where Qi~^^ is a heavy 6-quark field, v is its velocity, and Hj^^^ is the corresponding heavy 
meson field. This bosonization has to be compared with the matrix elements defining the 
meson decay constants fn {H = B, D) are the same when QCD corrections below mq are 
neglected (see 0,0]): 

= Jc-^) ' 

where Cy^^ are Wilson coefficients due to perturbative QCD for scales /i < mq {Q = b,c for 
H = B,D). We take /x = A^, which is the scale where perturbative QCD are matched to 
our hadronic matrix elements. 

For the VT-boson materializing to a D we obtain the bosonized current 

qll" Q^i^ — ^Tr l^h'-LHt^] , (17) 

where v is the velocity of the heavy c quark and ^ is the corresponding field for the D 
meson. 

For the b ^ c transition, we obtain the bosonized current 



Q\r>rLQ'J> -auj)Tr[m'^>j"LHl,-^>\ , (18) 

where C('^) is the Isgur-Wise function for the B ^ D - transition, and v' is the velocity 
of the heavy c-quark. Furthermore, uj = v-v' = v- v = Mb/{2M£)). Note also that from 
conservation of momentum we find the relation between the heavy quark velocities: 

Pb=Pd+Pd ^ = ^{v' + vy. (19) 

Mr 



Dt % 



FIG. 3: Non-factorizable chiral loops for — > DfD^ 



For the weak current for DD production (corresponding to the factorizable annihilation 
mechanism) we obtain 



(20) 



, -C(-A)Tr 

where \ = v ■ v' = [M|/(2M|,) — 1]. The Isgur-Wise function C(— -^) in (EOl) is a complex 
function, and not so well-known as for the b c transition. In the factorized limit, the 
matrix elements of the four quark operators are obtained by multiplying the bosonized 
currents above. 

In the following we will consider explicitely the decay mode 5° D~^D~ . The analysis 
of — s> D^D^ proceed the same way. To calculate the chiral loop amplitudes we need the 
(factorized) amplitudes for B'f —>■ D^D*^ and B^ —>■ D*^D*~ , which proceed through the 
spectator mechanism as in figure ^ The point is that the leading chiral coupling obtained 
from (|T2|l is a coupling between a pseudoscalar meson H, vector meson H* a light pseu- 
doscalar M (= 7c,K,r]). Using the bosonized currents in eqs. (fT7|l and p8|l . we obtain the 
following chiral loop amplitude for the process B^ — > DfDj from the figure El 



A{B^ 



where the factorized amplitude for the process -B" —>■ D^D^ is 



A{B^ ^ D+D-)f = fVcbV:, a2 C(^)/dMb v^M^ (A + u) 

v2 



(21) 



(22) 



The quantity R^ is a sum of contributions R12 from the left and right part of figure El 
respectively. In the MS scheme the results for i?^2 ^-re 



2(47r/)^ 



m 



2{Anf) 



9A 



(a; + A) ' 
0(^ + 1) 



Uln 



m 



K 



(A + cu) 



r(-A) - 1 hn 



m 



K 



A2 

X 



(23) 
(24) 
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FIG. 4: Non-factorizable contribution for DfD^ through the annihilation mechanism with 

additional soft gluon emision. The wavy lines represent soft gluons ending in vacuum to make 
gluon condensates. 



Adding these two contributions we find : 



m 



K 



{Anf) 



;9a 



[r{-uj) +r(-A)] - 1 Un 



A2 

X 



1 



(25) 



(UJ + X) 

As usual, the l/Nc suppression is due to p ~ Nc. The function r(x) is also appearing 
in loop calculations 0, 0] of the anomalous dimension in HQEFT (for a; > 1 and x < —1 
respectively) : 

1 



r{x) 



— - In (^x + Va;^ — 



r{—x) = —r{x) + 



ITT 



(26) 



which means that the amplitude gets an imaginary part. Numerically, we find 



0.12 - 0.26i 



(27) 



IV. NON-FACTORIZABLE SOFT GLUON EMISSION 

The genuine non-factorizable part (see eqs. (@]), (0) and Q ) can, within the framework 
presented in this section, be written in a quasi-factorized way in terms of matrix elements 
of coulored currents: 

{DiD;\Cwm%p = 8C2 {DtD-\-CLi,t-CL\G){G\dLin%Lm , (28) 

where a G in the bra-kets symbolizes emision of one gluon (from each current) as vizualized 
in figure |3 We observe that the annihilation mechanism amplitude in the non-factorizable 
case has the numerically favourable Wilson coefficient C2. 
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In order to calculate the matrix elements ()28|) . we will use a model which incorporates 
emision of soft gluons modelled by a gluon condensate. This will be the Heavy Light Chiral 
Quark Model (HLxQM) recently develo ped in This model belongs to a class of models 
extensively studied in the literature . For details we refer to ref. 

The Lagrangian for the HLxQM is 



0. 



HLxQM 



HQEFT 



xQM 



+ c 



Int 



(29) 



The first term is given in equation (fTT|) . The light quark sector is described by the Chiral 
Quark Model (xQM), having a standard QCD term and a term describing interactions 
between quarks and (Goldstone) mesons: 



xQM 



(30) 



Here m is the SU(3) invariant constituent light quark mass, and x is the flavour rotated 
quark fields given by xl = ^^11 , Xb. = ^Qr, where = {u, d, s) are the light quark fields. 
The left- and right-handed projections Ql and qn are transforming after SU{3)l and SU{3)ii 
respectively. In (jHUj) we have discarded terms involving the light current quark mass which is 
irrelevant in the present paper. The covariant derivative in ()30p contains the soft gluon 
field forming the gluon condensates. The gluon condensate contributions are calculated by 
Feyn,„an diagram techniques as in QMQQ. They may also be calculated by means 
of heat kernel techniques as in [15|, [161, [19 1 . 

The interaction between heavy meson fields and heavy quarks are described by the fol- 
lowing Lagrangian ^: 



C 



Int 



-G 



H 



X. 



(31) 



where Gh ~ V2m// is a coupling constant. In [4] it was shown how (fT^ could be obtained 
from the HLxQM. Performing this bosonization of the HLxQM, one encounters divergent 
loop integrals wh.ch wllin general be quadrafc. hnear- and logarithmic divergent Q. Also, 
as in the light sector jl^ the quadratic and logarithmic integrals are related to the quark 
condensate and the gluon condensate respectively. 

Within the model, one finds the following expression for the Isgur-Wise function 

2 



1+cu 



{1- p)+ pr{uj) 



(32) 



10 

where p is a hadronic parameter giving the deviation from the leading value P| : 

2m (1 + 3^^) + ^ 3 

Gh = ^P , P= , jVm^N ' f^ciH) = -mgiMH* - iW^n), (33) 

^ ^^^ + ^^> ^ 

where ri = (1 + 2/7r). Numerically, the deviation of p one is of order 10%. The simple 

expression in (jH^ is modified by perturbative QCD corrections down to fi = (analogous 

to eq. (HH)) and chiral loop corrections. Our model dependent expression ({u) in (jH^ give 

a good description of the Isgur-Wise function. 

The left part in figure |21 with gluon emision gives us the bosonized coloured current : 



(34) 



where is the octet gluon tensor, and -ff^^"* represents the heavy 5-meson fields. Similarly 
the (heavy) D- and Z)-mesons are represented by Hc^^ and ■* corresponding to a heavy 
quark field Q^J"^ and heavy anti-quak field Q^""* respectively. The symbol { , } denotes the 
anti-commutator. 

For the creation of a DD pair in the right part of figure El the analogue of (jH^ is 

Gl9s 



(35) 



X i-a^" + -{(T^^7■^} 

\n 4m(A - 1) ^ ' ' ^ 

where t = v' — v , and r = r(— A). Multiplying the currents in eqs. ()34p and ()35|) . and using 
the replacement: 

g'sG'^^Glf, Arr'{^G')^{g,^g.f, - g,pg,^) , (36) 

we obtain the bosonized version for the operator Qi in eq. © (see also eq. (j2HI)) as the 
product of two traces. (The expression may be simplified by using the Dirac algebra, but 
we do not enter these details here). 

Taking the pseudoscalar parts of (jH^ and we find the gluon condensate contribution 
for 5° — > DtD^ within our model: 



A(W DtD-,). ^ C^e^G') (l + if) . (37) 

(For our algebraic manipulations, the program FORM was useful). The ratio between 
this amplitude and the factorized one in (j^^ scales as Mo/iNcMs) times hadronic pa- 
rameters calculated within HLxQM. We define a quantity R'^ for the gluon condensate 
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amplitude analogously to in ()2ip and ()25|) for cliiral loops. Numerically, we find that the 
ratio between the two amplitudes in (jHTj) and is 

Rg ^ 0.055 + 0.16Z , (38) 

which is of order one third of the chiral loop contribution in eq. ()25|) . 



V. DISCUSSION AND RESULTS 



Our amplitude is complex as expected. In the chiral loop amplitude these are due to 
physical cuts (exchanges of physical particles) to the one-loop order we consider in this paper. 
The Wilson coefficients turn complex when the c-quark is treated within HQEFT. This 
is also the case for the matrix elements that these Wilson coefficients should be matched to. 
There is a potential problem with a quark model without confinement that the amplitude 
may get an imaginary part due to production of free quarks. Still, within HQEFT one can 
hardly distinguish from because of the reparametrization invariance. Thus, at the 
present stage, it is not clear how well our model describes imaginary matrix elements, and 
we will not go into such details here, as the numerical consequences turn out to be minor. 



Adding the amplitudes R^ and Rq and multiplying with the Wilson coefficient 
02 — 1.33 + 0.2z, we obtain the quantity: 



Rt = a2 (i?x + Rg) ^ 0.26 - O.lli . (39) 

Dropping the imaginary parts of the three quantities would give instead the value ^ 0.25. 
Anyway, we have found that the amplitude for is of order 15 — 20% of the 

factorizable amplitude for B'^ D^Dj, before the different KM-factors are taken into 
account. We obtain the branching ratios 

BRfm^DtD;) = 6.5 X 10-5 x|-^-^^^|2 (40) 
^ " ^ '0.041 0.223 0.25 0.9 ' ^ ^ 



and 



, y,, Rt C(^),2 



Si?(SO ^ D+D-) = 8.9 X 10"^ X 1—^ — ^ — ^ (41) 
^ ' ^ '0.041 0.974 0.25 0.9 ' ^ ^ 



The difference between the branching ratios is mainly due to the difference in KM factor. 
Taking into account the comments above, we end up with the conclusion that 

BR(B^ DjD;) ~ 7 X 10-^ , and BR{W^ DjD;) ~ 1 x 10"=^ . (42) 
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The ongoing searches at Belle might soon give the hmit on the rate S° — > DfDj, while 
the detection of the mode might be presently more difficult due to troubles with 
identification. 

S.F. thanks P. Krizan and B. Golob for fruitful discussion 
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